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Meson Green's functions and decay constants fp in different channels T are cal- 
culated using the Field Correlator Method. Both, spectrum and fp, appear to be 
expressed only through universal constants: the string tension a, a s , and the pole 
quark masses. For the S- wave states the calculated masses agree with the exper- 
imental numbers within ±5 MeV. For the D and D s mesons the values of fp(lS) 
are equal to 210(10) and 260(10) MeV, respectively, and their ratio f D Jf D =lM(3) 
agrees with recent CLEO experiment. The values fp(lS) = 182, 216, 438 MeV 
are obtained for the B, B s , and B c mesons with the ratio /b s //b=1-19(2) and 
/d//b = 1-14(2). The decay constants fp(2S) for the first radial excitations as well 
as the decay constants fy(lS) in the vector channel are also calculated. The dif- 
ference of about 20% between fp s and fo, fB s and fs directly follows from our 
analytical formulas. 



I. INTRODUCTION 



The decay constants fp in the pseudoscalar (P) channel, being important characteris- 
tics of mesons, in many cases can be directly measured in experiment, and therefore they 
can provide a precise manner to compare different theoretical approaches and check their 
accuracy. During the last decade the constants fp have been studied by many authors in 
potential models [1, 2, 3, 4, 5, 6, 7, 8], in the QCD sum rule method [9, 10], in lattice simula- 
tions [11, 12, 13, 14], as well as in experiment [15, 16, 17, 18, 19]. (The papers [2, 3, 4, 5, 6, 7] 
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contain references and a summary of numerous calculations of decay constants). 

The present article is devoted to the systematic derivation of the meson Green's functions 
in QCD and study of the decay constants for channels with arbitrary quantum numbers 
r, of which we specifically consider P and vector (V) channels. For the decay constant 
fr transparent analytical expressions will be obtained and in particular, using those, the 
difference between f Ds and fr>, fs 3 and can be easily explained. 

This paper is an improvement and extension of the earlier paper [20] devoted to the 
heavy-light (HL) pseudoscalars. Ref. [20] appeared before the systematic formulation of 
The Field Correlator Method (FCM) [21], in particular, before the derivation of the string 
Hamiltonian [22], therefore some steps in [20] were not accurately proved. In this paper 
we give a consistent and general treatment of the meson Green's function and its spectral 
properties. The main problem, which one encounters when addressing the spectral properties 
in QCD, is the necessity to include quantitative nonperturbative (NP) methods, which are 
responsible for the main dynamical phenomena: confinement and chiral symmetry breaking 
(CSB). 

In the FCM, introduced in [23], one derives the effective Hamiltonian, which comprises 
both confinement and relativistic effects, and contains only universal quantities: the string 
tension a, the strong coupling a s , and the current (pole) quark masses We use here 
the pole quark masses which correspond to the conventional current (Lagrangian) masses 
m q (m q ) [16]. The simple local form of this Hamiltonian, which will be called the string 
Hamiltonian (SH), occurs for objects with temporal scales larger than the vacuum gluon 
correlation length, T g 0.2 fm, i.e., it is applicable to all QCD bound states with an 
exception of toponium. Explicit calculations of masses and wave functions with the use of 
the SH have been done recently for light mesons [24] , heavy quarkonia [25] , and heavy-light 
mesons [26], and demonstrate good agreement with experimental masses, leptonic widths, 
and fine structure effects. 

As compared to QCD sum rules and lattice QCD this method has an essential advantage, 
because the radial and orbital excitations can be considered in this approach on the same 
grounds as the ground states. 

The calculation of spectral coefficients, like decay constants fr, needs an additional 
step, namely, besides using the SH the computation of all coefficients in the Green's func- 
tion, including the Dirac spinor structure etc. Moreover, for n and K mesons CSB is vitally 
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important. Recently the FCM was extended to include the effects of CSB [27], where it 
was shown that the phenomenon of CSB occurs due to confinement and two characteristic 
parameters of CSB - the chiral condensate and - were computed in terms of a. The 
calculation of /„- and fx can be done using a simple extension of general expression (23), 
derived in Ref. [27], while here we concentrate on the calculations of the masses and /r of 
HL mesons; our method also enables one to calculate fr for excited states and here decay 
constants will be calculated for the 2S states. 

One important technical problem, which is solved in this paper and allows one to calculate 
fr in all channels, is the accurate einbein reformulation of the Fock-Feynman-Schwinger 
Representation (FFSR), or the world-line representation, where the dynamical quark mass 
u q appears as an integration variable instead of the proper time. The previous step in this 
direction [22] has enabled one to write only the string Hamiltonian, but the whole Green's 
function was not attainable. Below in the FFSR we derive the explicit einbein form of the 
meson Green's function. 

The paper is organized as follows: in Section II the general einbein form of the meson 
Green's function is presented, while in Appendices A - C the details of the derivation are 
given. In Section III the masses of the heavy-light pseudoscalar mesons are considered. 
The decay constants are calculated in Section 4, while auxiliary variables are given in Ap- 
pendices D and E. Section V is devoted to the approximations used in our calculations and 
Section VI contains our concluding remarks. 



Here we derive the einbein form of the meson Green's function written as the path integral 
in the FFSR [28, 29, 30]. We start with the FFSR for the quark Green's function in the 
gluonic field A^, which contains both perturbative and NP contributions and Euclidean 
space-time is assumed everywhere: 



II. 



THE MESON GREEN'S FUNCTION IN THE FFSR 



S(x,y) = (m + D) 



-i 



(m - D)(m 2 - b 2 Y 



i 




(1) 
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with 



d 

D = D^, D^ = — - igAp, 
K = m*s + l - I (^Ur, 



i r fdz^ 

$<r(x,y) = Pexp (ig J A^dz^j exp J a^F^drj , (2) 



where 



(Dz) xy = hm M 



\a=b 



,ip(a-b) 



, Ne = s, 

T,k Az ( k )= x -y 



d A p 



(3) 



(2tt)' 

In Eq. (1) the role of the evolution parameter is played by the proper time s, whereas 
in (Dz) xy there is an integration over the fourth component z 4 (r), (0 < r < s), which is 
the Euclidean time of the particle. The crucial point now is to go over in Eq. (1) to the 
Euclidean time z^ = t as an evolution parameter. To get rid of the proper time s, one can 
use the so-called einbein method [31], which was applied to the FFSR in Ref. [22, 30], and 
here it is developed and used for the correlator of the currents. 

To this end the so-called dynamical mass (variable) uo(t) can be introduced via the relation 
between the proper time r and Euclidean time t: 

dT = ^ry G=!jT«(T)* = i|>(*). w 

The integrals in Eq. (1) can now be identically rewritten as 



d 3 A7(k) 



duj{k) 



(5) 



where 



( ) = ^^' l ^ k)= \~At ' NAt = x *-y* = T - ( 6 ) 



In Appendix A the representations (1), (5) are illustrated by calculating the free quark 
propagator, where the meaning of a; (A;) and lD in the momentum representation appears to 



be very simple: to = u(k) = a/p 2 + m 2 . 
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We now turn to the meson (quark-antiquark) case and consider the correlator Gr(x) of 
the currents jr(x): 

j r (x) = ^i(x)r^ 2 (a;), 

r = t a ® (1, 75, 7 M , 27^75) for S, P, V, and A channels, 

t a = y, tr(t a t b ) = U ab , (7) 

and 

G r (x) = (jr(x)j r (0)) v 

= AN C J Y r (D 3 z) x0 (D 3 z) x0 (Duj 1 )(Du 2 )exp(-K 1 - K 2 )W a . (8) 

In (8) we have defined the new quantity Y r : 

4F r = tr L {m x - D x )Y{m 2 - D 2 )F -> tr L {m x + uj x z)Y{m 2 - u 2 l)T, (9) 

which can also be written in the operator form: 

4F r = tr L {{ mi - ipi)r(m 2 + ip 2 )F), (10) 

with p$ (pi = p^fj) - the momentum of particle % as it is derived in Appendix B for the 
convenience of the reader (tri means the trace over Dirac indices). The resulting expressions 
for Y r in the V, A, S, and P channels are given below in Eqs. (35). 

In (8) the symbol W a = (x, y)& a2 (y, x) stands for the average value of the Wilson 
loop with the insertions of the operator 

A?> = exp (ga$ jT i^dr) = exp (go$ £ F^zfa))^^) (11) 

for the quark line and of the operator A?) for the antiquark line: 

= exp (-ga$ jT J^dr) = exp (-^g) jf * J^(*(f 2 ))_||-^ . (12) 

Since in HL mesons we will consider spin-effects as a perturbation, in the first approximation 
both factors A^ and A^ are replaced by 1 and W a is simplified. 

In Appendix C one can find all explicit steps for the derivation of the correlator (8) in 
the simplest case when the gluon interaction is absent, i.e., for W a = 1. One can see there 
that the quark-loop contribution is reconstructed with the correct coefficients. Now we turn 
to the case when the NP interaction is included in W a . 
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In general, W a contains all effects of the interaction which include: (i) the perturba- 
tive static gluon exchange; (ii) the radiative corrections to G r in the form of the operator 
anomalous dimension and corrections to /r; (hi) the NP contributions to Gr- To calculate 
all of them one can use the background perturbation theory [32, 33] and the FCM. For all 
hadrons of interest to us (with a size larger than T g 0.2 fm) the use of the FCM reduces to 
the appearance of the area-law factor in W a , which is accompanied by Coulomb, radiative, 
and spin-dependent factors. (For more discussion see the reviews [30, 34]). As was shown 
in Refs. [33, 34] this produces (apart from radiative corrections) the local interaction V(r), 
which for not so large angular momentum (L < 4) can be presented as W a = exp(— J Q T dtV), 
where the interaction, 

V(r) = V (r) + Vsd + AV string + Vse, (13) 
contains the static potential: 

V (r) = ar- , 14) 

3 r 

and the spin-dependent part Vsd(?") given by: 

Vsd = Vgs + V LS + V T , (15) 

with spin-spin, spin-orbit, and tensor terms; the self-energy contribution Vse [35], and also 
a "string correction" (occurring only for the states with L ^ [22, 24]). For the S'-wave 
mesons, considered here, 

V(r) = V (r) + V ss + V S E, (16) 

while for the spin-averaged masses only two terms are left in the potential, V(r) = Vq{t) + 
Vse- In Section V we shall also take into account radiative corrections and the operator 
anomalous dimension. 

One can now rewrite Gr in Eq. (8), separating cm. and relative distance coordinates, 

r) = Azi - Az 2 , p = ■ , uj t = ■ , 17) 

LV 1 +UJ 2 UJi + UJ 2 

where all coordinates are labelled with the index k as in Eq. (6). 

Integrating out the cm. coordinate dp and du+ = oj\ +u 2 , as shown in Appendix C, one 
arrives at the path integral in the relative coordinate drj, which can be expressed through the 
SH. Indeed, from the path integral formalism [30, 36] it is known that a general equivalence 
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relation holds: 



/(27T 



At/u r {k)f/ 2 



= (^|e- HT |y), 



where 



Taking into account the integral (see Appendix C): 



V^r(fc)' 

one obtains the important relation: 



(18) 



(19) 



(20) 



x = N r (0 



r e _HT 



(21) 



where the Hamiltonian is 



H = sjp 2 + m\ + ^p 2 + m 2 . + V(r). 



(22) 



Then with the use of the spectral expansion the expression (21) can be presented as 







r 6 -ht 



-M n T 



(Wl)„(w 2 )n 



(23) 



where y9 n and M n are the eigenfunction (e.f.) and the eigenvalue (e.v.) of the Hamiltonian 
(22). It follows from the extremum conditions (D.7) that the variables a>i and u;2 in Eq. (21) 
are defined as the operators, 

(24) 



p 2 + mf , 



while in Eq. (19) (u>j) n is the matrix element over this operator for the nL state. 

On the other hand for the l.h.s. of Eq. (21) one can also use the conventional spectral 
decomposition: 

d 3 P 



G r (x)d 3 x = I d 3 x(0\j r \n}(n\ 3r \0)e lP - 

n J 



x-M n T 



2M n (2ivy 



2M n 



Here we have used the standard definition for /p = fr(nS): 



(25) 



(0\j r \n,P = 0)=e r M n ft., 



(26) 
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where e r = for V and A channels, e r — 1 for S and P channels, while the polarization 



vector satisfies the normalization condition: 



E 4%H%) = ^- qj f- (27) 

fc=l,2,3 ^ 

Inserting the expression (23) into (21) and using the relation (25), one obtains the decay 
constant /f in the channel T for the nS state: 

2 = 2JV c (y r )|Pn(0) I 2 = 6|^(0)| 2 (Y T ) 

(0 1 } n (uj 2 ) n M n M n (wi) n (w 2 ) n " 

To derive the relations (23) and (28) we have used the essential property of the SH 

Eq. (22) that the average value (uJi(k)) of the operator Ui(k) is equal to the average of the 



operator a/p 2 + mf, and hence the average Ui = ^ ^2 k=1 Vi(k), denoted as (ui) n is 

(Qi) n = (ui(k)) n = (V^r^)n, (29) 
where the average is assumed to be taken over the eigenstate y2„(r): 



(uji) n = (ip n \y p 2 + mf\(p n ) . (30) 

In the most general case this average may differ from the path-integral average in 
Eqs. (4), (5), (8). However, the analysis of this problem, done in Ref. [37], shows that the 
difference between the two definitions is small (< 3% for the lowest states) and we assume 
here that the same accuracy holds for our basic relation (28). 

In Eq. (28) the mass of the 1 So- wave meson M„( 1 So) = M on — |Ahf + Ase, where M on 
is the e.v. of the Hamiltonian (22) with the static interaction (14): 



H = ^ml+p 2 + ^m 2 2 +p 2 + V (r), 
H (p n = M on ip n . (31) 

The general expression of the self-energy correction to the meson mass M(nS) is calculated 
in [35] (see Appendix E) (it comes from the NP contributions to the quark and antiquark 
mass) : 

ASE = £ { ) + Mui)[m + (ui)/2 + e(u)]*} ' (32) 

where the factor r]f(i) depends on the flavor of the i-th quark (antiquark) and its analytical 
expressions (E.2) are deduced in Ref. [35]. For the u(d), s, c, and b quarks the values 

7] u{d) = 1.0, r] s = 0.65, r] c = 0.35, r] b = 0.03 (33) 
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are obtained in Appendix E 1 . 

It is worthwhile to notice that for the b quark i] b is small and its contribution to A SE is 
small (~ 1 MeV) and can be neglected. Therefore, for the B, B S) and B c mesons we have 
to use in Eq. (32) only a contribution which comes from the lighter quark (antiquark) q\{qi) 
denoted later by the index 1. For the D and D s mesons both terms (with % — 1 and i = 2) 
are taken into account, although the contribution from the c-quark term is small, around 
-20 MeV. 

Thus, the scheme of our calculations of / r is as follows: 

1. First, one calculates the e.v. M on of the SH (22) with the static interaction (14) and 
then takes into account the self-energy (32) and the HF corrections: 

M n ( 3 Si) = M on + A SE + ^A HF , M n CSo) = M on + A SE - ^A HF . (34) 

2. The values (u>i) and (u> 2 ) in Eqs (28) and (32), are the matrix elements (m.e.) of the 
kinetic energy term defined in Eq. (24). 

3. The factor (Yr) (in the channel T) can be computed in terms of the momenta of a 
quark and an antiquark, or in the cm. system in terms of the relative momentum p, with 
the following results for the m.e. (Yr) (see Appendix B): 

(y v ) = mi m 2 + (u)i){u)2) + g(P 2 )> 

(Yp) = mi m 2 + (tfiXwa) - (p 2 ) = (Y A4 ). (35) 
Here we use for the operators Yr the notations: 

Y v = ^^tr[(m! - D^iirrii - D 2 hi\, 

i 

Y P = Y Aa = -[trim! - ^1)7475(^2 - ^2)7475]- ( 36 ) 

In the case of the P channel with both mi,m 2 — > due to CSB there appears an additional 
mass term in Eq. (35), which can be computed through field correlators. (The ir and K 
mesons will be considered later [27]). 

For the calculations of (Yr) one needs also to know the m.e. (p 2 ) n ,s and the w.f. at the 
origin, (p n (0) = -R„(0)/ \ZAtt; they are given in Appendix D. 

1 Note, that these values of r]f(i) i n (33) are not fitting parameters, but calculated through the same input 
rrii and a 
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4. Our calculations are done with the relativistic SH (22) which was derived in einbein 
approximation (EA), therefore the w.f. at the origin must also be calculated with the use of 
the EA (see Appendix D), which provides an accuracy of < 5% [26, 37]. In the nonrelativistic 
limit, vrti ^> ^fa, one can easily find that (pi) ~ mj, while (p 2 ) ~ 0(a) can be neglected in 
Eqs. (35), and therefore 

(F v )nr ~ 2m im2 + 0(a), (F 5 )nr ~ 0(a), 
(Ya^nr ~ 0(a), (Ya 4 )nr = 2mim 2 + 0(a), 

(Yp) nr = 2 mi m 2 + 0(a). (37) 

Thus in the nonrelativistic limit for in the V and P channels one obtains the well-known 
result [1]: 

(f?) 2 NK = ^\M0)\ 2 (T = V,P), (38) 



n 



while in the S channel fs — > 0. 

5. As a final step one needs to compute the radiative corrections to fp, which come 
from the short-distance (large momentum) perturbative gluon contributions. Neglecting 
interference terms they can be written as in [20, 30], 

(W <r > = (WbGE)(W non pe rt > (39) 

with 



/w \ ry I 4 f [ dz±dz' 4 a s (z - z') \ 

(W OGE ) = Z m exp ( -— / / — — ) , (40) 



where Z m is a regularization factor. After separating the Coulomb interaction in H in this 
way, one gets the correction to (W a ), and f 2 can be written in the form: 

ft^ivfl iv=l + c v a s + 0(al). (41) 

Another important contribution from perturbative gluon exchanges (GE) is the account of 
Asymptotic Freedom (AF) in the coupling constant a s in (40), which is especially important 
for the value of (p n {0) in the 5- wave channels. In our calculations we use the GE interaction 
where in the strong coupling a_e(r) the AF behavior is taken into account. Nevertheless, it is 
of interest to compare the w.f. at the origin with and without AF behavior in the GE term, 

2 

introducing the factor paf = < ^n AF ^(0)/v? re (0) , which appears to be around 0.80 for the IS 
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TABLE I: The pole quark masses m q (pole) (in GeV), used in this paper, and the constituent 
masses m q from [2, 4, 5]. 



quark 


b 


c 


s 


u 


d 


m q (pole) this paper 


4.78 


1.40 


0.180 


0.005 


0.008 


m q from [4] 


4.655 


1.511 


0.216 


0.071 


0.071 


rriq from [2] 


4.977 


1.628 


0.419 


0.220 


0.220 


rriq from [5] 


5.158 


1.755 


0.535 


0.371 


0.377 



states [20]. Then /p (with AF taken into account) can be expressed through /f? (where AF 
is neglected) as /j 2 . = /r frPAF- 

6. We conclude this section with the discussion of the input parameters rrii, a s and a. 
We take a — 0.18 GeV 2 for all HL mesons (as in light mesons and in heavy quarkonia [38]); 
mi are the conventional pole masses which are defined through the Lagrangian (current) 
masses in the MS'-scheme (see [16, 39] and references therein): 

rm = ™mK™ms) {l + l^Mll + ^ + 0(Qb) s J . (42 ) 

We use here the pole masses presented in Table I. They correspond to the conventional 
current masses m c = 1.18 GeV and m b = 4.20 GeV while for the strange quark the pole 
mass, m s = 180 MeV is taken at the scale /x p=j 0.5 GeV (r ps 0.5 fm) and therefore it should 
be larger than the standard m s (2 GeV) ps 100 MeV [16]. For the u and d quarks m u = 5 
MeV, m d = 8 MeV are taken. 

One can see the essential difference between the input masses in our SH and the con- 
stituent masses used in the relativistic instantaneous Bethe-Salpeter Method [5], and also 
in the spinless Salpeter equation [2]. The b, c, and s quark masses from [4], where the 
relativistic Dirac Hamiltonian is used for a light quark, are not large and can be considered 
as the pole masses, while m u = m d = 71 MeV seem to be too large. 

7. The coupling a st (r) in the GE term is taken here as the vector coupling as(r) in 
background perturbation theory (in two- loop approximation) from [39]. At large distances 
this theory predicts saturation: as(r) — > a cri t = const (this coupling as(r) is in good 
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agreement with lattice potential at small distances [40]): 

2 f°° singr 

a B {r) = - dq a B {q), 

n Jo Q 

^^^i 1 -^/' * fl = ln T|- (43) 

Here the background mass Mb = 1 GeV can be expressed through the \fa [38], while the 
QCD constant is given by 

31 _ 10 „ 
9 

2 A) 



A B = AM5exp(^-/^). (44) 



It is of interest to notice that in our calculations of HL meson properties the value rif = 3 
is strongly preferable. 

In this way the problem is uniquely defined and no fitting parameters are introduced. The 
resulting masses and decay constants of HL mesons appear to be unbiased theoretical pre- 
dictions which will be compared with recent experiments, lattice data, and other theoretical 
predictions. 

III. MASSES OF HEAVY-LIGHT MESONS 

In the relativistic SH a correct choice of the static interaction Vo(r), which defines the e.v. 
of the unperturbed Hamiltonian H (31), is of great importance. Here, for all HL mesons 
we take the static potential from [38, 39]: 

V B {r)=ar-\^- (45) 
3 r 

with a = 0.18 GeV 2 , the number of flavors n/ = 3, the vector coupling as(r) (43) also 
contains 

A B (n f = 3) = 360 MeV, M B = 1.0 GeV. (46) 

For this choice the saturated (critical) value of the vector coupling, reached at large r, is 
equal a crit (n/ = 3) = 0.495, while at r ps r = 0.5 fm, a B (r ) = 0.43. 

It is very convenient to start the calculations with the spin-averaged masses M cog (nL), 

M cog (nS) = M (nS) + A SE (nS), 

M cog (nL) = M (nL) + A SE (nL) + A str (nL) CM 0), (47) 
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TABLE II: The spin-averaged masses M cog (lS) (in MeV) a) . 



Meson multiplet 


M (IS) 


A SE 


M cog (15) 


M cog (exp) 6 ) 


D-D* 


2139 


-164 


1975 


1974.8(4) 


D s ~ D* s 


2177 


-105 


2072 


2076.0(6) 


B-B* 


5433 


-120 


5313 


5313.5(6) 


B s ~ B* s 


5468 


-72 


5396 


5400.7(32) 


Bc-B* 


6332 


-17 


6315 




a ) The Ase are calculatec 


in Appendix E. 



6 ) The experimental numbers are taken from PDG [16] and [41]. 

since these masses do not depend on any additional parameter. In contrast to M cog the 
masses of the singlet and triplet states depend also on the parameters defining the HF inter- 
action, or for states with L ^ the masses M(J PC ) depend also on the coupling «fs(/^fs)) 
which defines the fine-structure splittings, and these couplings are not well determined. In 
Eq. (47)M (nL) is the e.v. of Eq. (31) (see Tables II,V).The self-energy term A SE (nL) (32) 
and the string correction A str (nL) = (H stT ) (D.12) are defined by analytical expressions and 
discussed in Appendices D and E. 

The calculated M cog (lS), as seen from Table II, agree with the experimental num- 
bers with an accuracy better than 5 MeV. It is of interest to notice that the difference 
M cog (Df) — M cog (D + ), which in experiment is 99 MeV, is only partly due to dynamical 
reasons: M (D S ) — M (D) = 38 MeV, but mostly occurs due to the difference in the SE 
contributions, equal to 60 MeV. 

With the use of our number M cog (B c ) = 6315 MeV and M(B c ) cxp = 6275(7) MeV [41] 
one can predict the mass of the vector state B*(l 3 Si): 

M(B* C ) = 6328(7) MeV, if M(B C ) = 6275(7) MeV. (48) 

The masses of the triplet and singlet HL mesons are calculated taking into account the 
HF interaction which in general contains both perturbative and NP contributions: 

A HF (n 1 S) = AP F + A^, (49) 
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where the P term with the one-loop correction is taken from [42]: 




Here p = ^/3q — § — f In 2 and for n/ = 3, p(nf = 3) = 0.5635, so that the one-loop correction 



is about 6%. From here on for simplicity we use the notation u>i,u>2 instead of (o>i) n , (&2)n- 
The NP contribution to the HF splitting was derived in [43] and with good accuracy it 
is given by the m.e.: 



Here G 2 is the gluonic condensate, for which we take the value G 2 = 0.043 GeV 4 [43], which 
provides the correct value of the string tension a = 0.18 GeV 2 . For the D and D s mesons the 
NP term (51) appears to be not small (as well as for the J/tp — r/c splitting): around 10 MeV 
if the gluonic (vacuum) correlation length is T g = 0.2 fm [44]. However, at present the value 
of T g is not known with high accuracy and T g = 0.28 fm was obtained from unquenched 
lattice data [45]. Due to this uncertainty in T g we have two possibilities to describe the HF 
effects in HL mesons: 

A. In the case of T g m 0.2 fm, A^p ps 10 MeV is not large. Then in the perturbative term, 
Eq. (50), we need to take a rather large coupling: a H F = 0.40 for the D and D s mesons and 
q;hf = 0.32 for the B and B s mesons, to reach agreement with experiment. 

B. In the case of large gluonic length, T g = 0.3 fm, the NP contribution A^p appears 
to be larger, about 22 MeV, and therefore the coupling q;hf can be taken smaller. In this 
case for the D and D s mesons the value «hf = 0.365 (rather close to that for the J/\I/ — r] c 
splitting [46]) gives rise to agreement with the experimental HF splitting. In Table III the 
calculated HF splittings both for T g = 0.2 fm and T g = 0.3 fm are given. 

Taking the HF splittings from Table 3 the masses of the singlet and triplet states can be 
calculated; they are presented in Table IV (M cog (15') are taken from Table II). 

The calculated triplet and singlet masses for the ground states appear to be in good 
agreement with the experimental numbers (with an accuracy < 5 MeV). In our analysis we 
have observed that for the D and D s mesons the value M cog (lS) is very sensitive to the pole 
mass of the c quark and high accuracy can be reached only for m c (pole)= 1.40 GeV; at the 
same time we cannot exclude that for the B c meson the choice of m c = 1.38 GeV is also 
possible. 




(51) 
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TABLE III: The hyperfine splittings (in MeV) in HL mesons. Case A: T g = 0.20 fm, a H F = 0.40 
for the D and D s mesons and «hf = 0.32 for the B and B s mesons. Case B: T g = 0.30 fm, 
ckhf = 0.365 for the D and D s mesons and «hf = 0.305 for the B and B s mesons. 





T g = 0.20 fm 


Multiplet 




aNP 
^HF 


A HF (tot) 


A H f (exp) 


D* - D 


131.6 


8.9 


140.5 


140.64±0.10 


D* s - D s 


130.7 


8.9 


139.6 


143.8±0.4 


B* - B 


43.5 


1.5 


45.0 


45.78 ±0.35 


B* s - B s 


44.2 


1.4 


45.8 


44.2±1.8 




T g = 0.30 fm 


D* - D 


119.4 


21.8 


141.2 


140.6±0.1 


D* s - D s 


118.6 


21.6 


140.2 


143.8±0.4 


B* - B 


41.3 


3.7 


45.0 


45.8 ± 0.4 


B s -B 


42.0 


3.4 


45.4 


44.2±1.8 



TABLE IV: The masses (in MeV) of the ground states (for T g = 0.3 fm). 





D± 


D *± 


D s 


D* 

^ s 


B 


B* 


B s 


B* s 


this paper 


1869.1 


2010.3 


1966.8 


2107.1 


5279.3 


5324.2 


5362.0 


5407.4 


Expert 


1869.3 
±0.4 


2010.0 
±0.4 


1968.2 
±0.4 


2112.0 
±0.6 


5279.0 
±0.5 


5325.0 
±0.6 


5367.7 
±1.8 


5411.7 
±3.2 



a ) The experimental numbers are taken from PDG [16] and M(B*) from Ref. [19]. 

For the B c meson, with the use of q;hf = 0.26 and M cog = 6313 MeV from Table II, we 
obtain A RF (B C ) = 44 MeV, which gives M(B C ) = 6280 MeV close to M{B c )(exp) = 6275(7) 
MeV [42] and M(B*) = 6324 MeV. 

Finally, in Table V we give the calculated masses M(2S) for the first radial excitations, 
which are not yet found in experiment. Our prediction for the singlet and triplet masses 
of the radially excited HL mesons strongly depends on the value of «hf taken. If for the 
D(2S) and D S (2S) mesons one takes the value ojhf^) = 0.30, as in the case of r) c (2S) [46], 
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TABLE V: The masses M cog (2S), M(2 1 S ) and M(2 3 S 1 ) (in MeV) for heavy-light mesons (a H F = 
0.30). 







D{2S) 


D S (2S) 


B(2S) 


B S (2S) 


B c {2S) a ^ 


this paper 


M (25) 


2758 


2797 


5998 


6034 


6868 




M cog (2S) 


2615 


2702 


5888 


5966 


6852 




M(2 1 S ) 


2560 


2646 


5864 


5941 


6821 




M{2 3 S l ) 


2633 


2721 


5896 


5974 


6862 


from Ref. [4] 


M(2 1 S ) 


2589 


2700 


5886 


5985 






M(2 3 Si) 


2692 


2806 


5920 


6019 




from Ref. [2] 


M(2 1 S ) 


2580 


2670 


5900 


5980 


6855 




M{2 3 Si) 


2640 


2730 


5930 


6010 


6887 



a) For the B C (2S) meson we use the value «hf = 0.26. 



then the values M(2 1 Sq) and M(2 3 Si), given in Table V, are obtained. 

The triplet and singlet masses, calculated here, are rather close to those from [2, 4], 
nevertheless for the 2 1 S states our numbers are systematically lower by ~ 30 — 50 MeV. 

In this paper we do not consider orbital excitations of HL mesons, it will be done in our 
next paper. Still, we would like to notice that our values of M(B sl ) and M(B s2 ) lie in the 
region 5.82 — 5.83, i.e., approximately by 100 MeV lower than in Ref. [4] and close to the 
recent experimental data [19]. 



IV. DECAY CONSTANTS OF HEAVY-LIGHT MESONS 

The general formula for the decay constants /p (23) can be rewritten (later on we shall 
use the notation u>i,u>2 for (u>i), (u 2 ) ) as follows 

= iSti* < 52) 



which contains the w.f. at the origin R n (0) = y? n (0)/v47r, the average kinetic energies 
cji,cj 2 , the meson mass M n = M(nS), and also the m.e. (p 2 ) in the factor (Yp). All these 
auxiliary values are given in Appendix D (Tables IX and X). Then the values of f P can be 
easily calculated. They are given in Table VI, together with the experimental, unquenched 
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TABLE VI: The decay constants / P (in MeV). 





D 


D 8 


B 


B s 




Ref. [20] 


206 


252 


174 


— 


- 


Ref. [7] 


230(25) 


248(27) 


196(29) 


216(32) 


322(42) 


Ref. [6] 


234 


268 


189 


218 


- 


lattice [11, 14] 


235(22) 


266(28) 




206(10) 




quenched 












lattice [11, 13] 


201(20) 


249(19) 


216(38) 


259(32) 


440(2) 


n/ = 2 + 1 












this paper 


210(10) 


260(10) 


182(8) 


216(8) 


438(10) 


experiment 


222.6(20)") 


280(23)°) 


160^0 b) 
229(70) c ) 







a ) The experimental values are taken from Ref. [17 . 
6 ) BaBar data [18] 
c ) Belle data [18] 



and quenched lattice data, and some other theoretical analyses. 

From Table VI one can see that our central value of fs d is 15% smaller than the one in 
unquenched lattice QCD [13], but rather close to /# in relativistic models [6, 7]. 

For the analysis of experimental data on direct measurements of the leptonic decay, 
P — > lis, it is important to know the ratios of the decay constants, which in our calculations 
are, 

^ = 1.24(3), ^ = 1.19(3), ^ = 1.20(3); ^ = 1.15(2). (53) 
Id Jb jb s jb 

These ratios are in good agreement with recent lattice data (unquenched) [13, 15] and close 
to the experimental number obtained by the CLEO collaboration with f D j ' f D = 1.27(14) 
[17]. It is of interest to compare these ratios with other theoretical calculations, which are 
typically smaller than our numbers, and also with recent lattica data (see Table VII). 

Due to the large theoretical errors in the ratios (d = /d s //d an d Cb — IbJ Jb (see 
Table VII) one cannot judge what is the true value of Cd and (b, however, from our general 
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TABLE VII: The ratios f D Jf D , f B Jf B , and f D Jf B .. 





JdJJd 


IbJ/b 


IDs/ fB^ 


RPM [6] 


1.15 


1.15 


1.23 


BS [7] 


1.08(1) 


1.10(1) 


1.15(1) 


lattice [13] 
unquenched 


1.24(8) 


1.20(4) 


1.01(8) 


this work 


1.24(3) 


1.19(3) 


1.20(3) 


experiment 


1.27(14) 







a ) The ratio of the central values is taken for fu a and f Bg ■ 

formula (52) for f P it follows that f Ds and f D as well as f Bs and f B have to differ by ~ 20%. 
This happens due to the presence in Eq. (52), through (Y p ), of a term proportional to mim 2 : 
mi = 0.008 GeV for D + and m\ = 0.18 GeV for Df, m 2 = 1.40 in both cases. If one neglects 
this term then the value fo d = jo u = 208 MeV practically does not change, while for the D s 
meson one obtains the essentially smaller number 215 MeV, instead of fp s = 260 MeV in 
the case with mi = m s = 0.18 GeV. Thus, the ratios (d and ( B appear to be very sensitive 
to what one takes for the pole mass of the s quark and can be used as a convenient criterium 
to choose m s at low renormalization scale. Note that the factor |-Ri(0) | 2 / (u; 1 c^ 2 ) = 0.35(1) 
turns out to be the same both for the D and D s mesons. Our number for ( D = 1.24(3) (with 
m s = 0.18 GeV) is in agreement with the experimental number £o(exp) = 1.27(14) [17]. 

The first radial excitations of HL mesons are considered here neglecting open channels 
which can decrease the w.f. at the origin. In this approximation 

f P (D(2S)) = 167 MeV, f P (D s (2S)) = 201 MeV, 

fp(B(2S)) = 168 MeV, f P (B s (2S)) = 194 MeV, f P (B c (2S)) = 347 MeV. (54) 

From these numbers it follows that for the 2S states /d s //d — 1-20 and f B Jf B = 1.15 are 
changed only by about ~ 4%, although the absolute values of f P (2S) appear to be ss 20% 
smaller than for the ground states. 

One can compare our results for f P with other theoretical calculations (see Table VI). The 
agreement with potential-model results [6, 7] is evident with the only exception for f P (B c ): 
our number f P (B c ) = 439 MeV is about 30% higher, and is close to the one calculated in 



19 



TABLE VIII: The vector decay constants fy (in MeV) and the ratio /v//p for the IS* states in 
heavy- light mesons. 





D* 


D* 


B* 


B* 


K 


(p 2 ) in GeV 2 


0.273 


0.291 


0.359 


0.383 


0.784 


fv [5] 


340(23) 


375(24) 


238(18) 


272(20) 


418(24) 


fy, this work 


273(13) 


307(18) 


200(10) 


230(12) 


453(20) 


fv/fp [5] 


1.48(26) 


1.51(26) 


1.21(27) 


1.26(28) 


1.30(24) 


fy/fp this paper 


1.27(5) 


1.17(4) 


1.08(4) 


1.07(4) 


1.03(3) 



unquenched lattice QCD, where fp(B c ) (lattice) = 420 (20) MeV in [13]. A detailed analysis 
of fs c (in the framework of the potential model and QCD sum rule approach with radiative 
corrections taken into account) [47] gives the value /b c ~ 400 MeV, which is approximately 
10% lower than our number. 

Finally, some remarks about decay constants in V channels. As seen from the expressions 
for (Yy) and (Yp), fy has to be larger than fp (since My and M P are close to each other), 
while (p 2 ) enters with different signs 

f 2 (nS) = [m x m 2 + uj^ 2 + \(p 2 ) nS ) M P (nS) 

f 2 (nS) (mim 2 + WlW2 -(p 2 ) nS ) M y (nS) ' 1 j 

These ratios (for the ground states) are given in Table VIII. For the heavy B c meson this 
ratio is approaching unity. 

From Table VIII one can see that our numbers for the vector decay constants are system- 
atically lower than fy (central values) from Ref. [5]: by ss 20% for the D and D s mesons 
and by 15% for the B and B s mesons, although their values lie within the large theoretical 
errors. 



V. APPROXIMATIONS 



We discuss here the approximations we made and the accuracy of our results. The starting 
expression for the current correlator (8) is exact, because the FFSR is an exact representation 
of the meson Green's function. The main approximation refers to the transition from the 
FFSR path-integral to the local Hamiltonian formalism (19) neglecting quark pair and hybrid 
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production and taking spin-dependent interactions as a perturbation. The accuracy of this 
approximation is determined by several factors: 

1) Neglect of hybrid excitations in the Wilson loop, which actually leads to a multichannel 
Hamiltonian. The hybrid admixture was shown to be small, of the order of few percent for 
ground state mesons [48, 49]. 

2) The use of the EA to define the w.f. at the origin and some m.e. The corresponding 
accuracy was checked in [37] and shown to be around 5%. 

One should stress that the SH contains only well defined fundamental parameters: the 
pole quark masses, Aqcd, and the universal string tension a, and does not contain any fitting 
parameters, in particular, there is no overall constant often used in the static potential, or 
in the mass. The accuracy of the SH was checked for light mesons [24], heavy quarkonia 
[25], heavy-light mesons [26], glueballs [50], and hybrids [49]. In all cases meson masses are 
in agreement with experimental and lattice data with an accuracy of a few percent. 

Concerning the decay constants fp, a special sensitivity occurs in the w.f. at the origin 
|y? n (0)| 2 , i.e., to the behaviour of the GE potential (or the vector coupling a st (r)) at small 
distances, not only in the AF region (r < 0.1 fm) but also in the region 0.1 fm < r < 0.3 
fm, which in [39] called the intermediate region. This behavior is known quite well for the 
perturbative part, where Aq CD is well known (for rif — 5 A^ = 217^23 MeV [16]). The 
major uncertainty comes from two sources: 

(i) the behaviour of the spin-dependent part, in particular, the HF interaction, where 
smearing of the 5 function can drastically change the wave function for systems of small 
size, R < 0.4 fm. For HL mesons with R > 0.6 fm this effect is becoming less important 
(for more discussion of the influence of the HF interaction see [51]). 

(ii) We estimate the accuracy of the resulting |(/? n (0)| 2 better than ±5%. Thus the accu- 
racy of our computed values of fp is expected to be < 8%, while the ratio of decay constant 
has better accuracy, < 4%. 

At this point one should discuss the effects which were unaccounted for till now. First 
of all, this concerns the radiative corrections due to transverse gluon exchanges and higher 
loops. These corrections contain the (pseudo) evolution factor, governed by the operator 
anomalous dimension, and considered in [52] two decades ago, 




(56) 
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One can estimate that when going from B to D mesons this factor changes only by 4%. 

Nowadays the radiative corrections are done within HQET and known to three-loop 
accuracy [53], e.g. for fo one has 

,(3) [ o$ (m c ) 



1 



1 + Cr 



(4)/ 



4% 



+ c\ 



4:71 



+ 



F®(m e ) + O ( , (57) 



where F^(m c ) is subject to the operator anomalous dimension correlations, and C\ 
etc. As a result, e.g. for the ratio one has [53] 



Id V m b 



1 + rAx 



1) ^ 4) K) +Q 



4tt 



(4)' 
47T 



+ 4 



A-QCD ^ 



-2C f 



(58) 



with r x = ||C2 + ff^f = 3.58, and x w 1.56. One can see that loop corrections contribute 
less than 4% and can be neglected within the accuracy of HQET and our approximations. 

From Eq. (58) it follows that in HQET the ratio /b//b < 1 as it happens in relativistic 
models [2, 5, 6, 7] and also in our calculations where /b/ Id ~ 0.87(2), while in unquenched 
lattice data this ratio is larger than unity, but has large computational error (~ 20%). 

We now turn again to Table VI and discuss our results in comparison to other calcula- 
tions. The first important point, which should be stressed, is that our input is minimal and 
fundamental, e.g. the pole masses in Table I correspond to the current masses quoted by 
PDG [16]. This is in contrast to many relativistic potential models with spinless Salpeter [2] 
or Bethe-Salpeter Hamiltonian [5, 7], where the constituent quark masses are used as input. 
In addition an overall constant is usually introduced in the interaction. In many respects 
our approach can be considered on the same grounds as lattice simulations, or the QCD- 
sum-rule approach (where instead of a several additional condensates are used). Moreover, 
the advantage of our approach is that excited states can be considered as well as the ground 
states. For excited states a quark (antiquark) of a given nL state has its characteristic 
"constituent" mass {uj)„,l, which grows for higher nL. Due to this effect, and also to the 
negative string correction, in our approach the masses of the radial and orbital excitations 
appear to be smaller (~ 20 — 40 MeV) than in other relativistic models (see Table V). 

Looking at Table VI one notices a good agreement (within 6 — 10%) of our results with 
unquenched lattice data and experiment (for fo and /dJ. Thus one can conclude that the 
Field Correlator Method and its essential part, the effective String Hamiltonian, appears to 
be successful in the prediction of / r as well as in other tests done so far [21, 22, 23, 24, 25, 26]. 



22 



VI. CONCLUSIONS 

In this paper we have studied the current correlator Gr(x) = (jr(^)ir(O)) and the in- 
tegral Jr — J Gr(x)dx in an arbitrary channel T with the use of the FFS path-integral 
representation. This method allows one to express decay constants fr(nS) for HL mesons 
through well-defined characteristics of the relativistic SH which was successfully used before 
in light mesons and heavy quarkonia. 

It is essential that the SH does not contain any fitting parameters, being fully defined by 
universal fundamental values: the conventional pole masses, the string tension a, and the 
strong vector coupling a st (r). 

The analytical expressions, obtained here for / P , show that the decay constants fo s and 
f Bs for the ground states have to be always larger by 20-25% than f D and f B because of 
the large difference between the pole (current) masses of the strange and light u(d) quarks. 
This theoretical statement is supported by recent experimental data [17]. 

In our analysis we have observed that 

• The calculated masses M(l 1 5'o) and M(l 3 Si) of all HL mesons agree with the ex- 
perimental numbers within ±5 MeV. Our prediction for B* is M(B*) = 6325(10) 
MeV. 

• The calculated masses of the first radial excitations M{2 1 Sq) for the D and D s mesons 
appear to be ~ 40 MeV lower and for the B, B s , and B c mesons ~ 20 — 30 MeV lower 
than the numbers from Refs. [2, 4]. 

• For the decay constants the values f D = 210(10) MeV and f Ds = 260(10) MeV are 
obtained. Their ratio /d s //d — 1-24(3) is close to the experimental number 1.27(14) 
[17]. 

• Our decay constants f B = 182(8) MeV and f Bs = 216(8) MeV give the ratio f B J f B = 
1.19(3) which agrees with a recent unquenched lattice number [11, 13]. 

• In the V channel the ratio fv/ fp is monotonically decreasing while going from the D 
meson to the heavier mesons: it is equal to 1.27(6), 1.17(4), 1.08(4), 1.07(4), 1.03(3) for 
the D, D s , B, B s , and B c mesons, respectively. For the D* and B* mesons calculated 
here, fy turned out to be ~ 20% smaller than in Ref. [5]. 



23 



APPENDIX A: FREE QUARK PROPAGATOR IN THE EINBEIN 
PATH-INTEGRAL REPRESENTATION 

One starts with the FFSR for the free quark propagator, which can be written as 

POO 

S(x,y) = (m-d) ds(Dz) xy exp(-K), (A.l) 
Jo 

and introduces the einbein variable, or the dynamical mass, u(t) as in (4), such that the 
function K can be rewritten as 

K = m s + - V ' dr 



4 Jo V dr 

In (Dz) xy , Eq. (3), there is an integration over the time components of the path, namely, 

< D *)-n^§'(E A *- 1 '). (a.3) 

k 

where T = x 4 — y 4 . With the use of (4) one can rewrite the integration element in (A.3) as 
follows (t = z 4 ), 



dAzAk) rtI .,. fe du(k)VAt n \ At . , . 



Moreover, the (^-function in (A.3) acquires the form 

6 (J2 Az4 ~ T ) =S(2us-T), (A.5) 

where we have defined 

i r 

iu = - Lu( T )dT. (A.6) 
s Jo 

As a result in (A.l) one can integrate over ds using the (^-function (A.5), and rewrite ds(Dz) xy 
as it is shown in Eq. (5) of the main text. 

Then one can write the Green's function as follows: 

S(x,y) = (m-d)JH g^ e -* ^*> ( ^ e^^», (A.7) 

where K is given in (A. 2), and /, in Eq. (6). 
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The integration over d 3 Azi(k) yields 

S(x,y) = (m-d) J j—^e V (A.8) 
Taking into account the relation, 

r°° Mki e ^ (Wh4^) = [E e -**V*^ , A 9 x 

one has for the scalar part G{x,y), defined by S = (m — d)G, the following expression: 

r ^3 pip-(n-v)-fo <&\J P 2 +m 2 

G(x,y) = / n > ( A -!0) 



m 2 



(2tt)3 2v ^T 

where we have used the relation following from the stationary point in the integral (A. 9): 

N 



X - jf co(r)dr = ^ = vV + m 2 . (A.ll) 



The expression (A. 10) can be compared with the conventional integral, 

j^U-o 5 ~, r = x-y, (A.12) 

which reduces to (A. 10) after integrating over dp 4 for T > 0. The expression (A.12) is the 
standard form of the free propagator. 



APPENDIX B: CALCULATION OF THE FACTOR (Y r ) 

Here we use and extend Appendix 1 of Ref. [20], considering the quark propagator (1) 
in the gluonic field A^. S(x,y) = [m — 7 M (d/dx^ — igA^)] G(x,y). The derivative d/dx^, 
acting on G(x,y), differentiates only the 5-function, 
5 (x^ — y^ — ^2k=i an( i can ^ e rewritten as a derivative in Az^(n): 

{'-«-£ Az(k) ) = -BE^m 1 { x - y -t Az(k) ) ■ (R1) 

Integrating by parts in the expression for G(x,y), one obtains 

D ll G(x, y) = J™ dse- ms (Dz) xy e- K $ a (x, y) (-^^ + 0(y/i)^ . (B.2) 
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In the limit e — * 0, iV — > oo one has A ^ fc ) — > and using the relation (4), one 



2e 2dr 

obtains 



T = S 



dz (t) 

D,G{x,y) = -u(r = 8 )-£±G(x,y), (B.3) 

where the r.h.s. of (B.3) is a symbolic writing implying that oj^- should be under the 
integral in G(x,y). Finally, when G(x,y) is expressed via the Hamiltonian H u (as in the 
definitions (19) and (21)), then one realizes that UiZp\t) = p$ and finds the relation: 

(m - D)G = (m - ip)G, (B.4) 

which will be used throughout the paper. 

As a check one can see that (B.4) yields the correct form of the free quark propagator in 
Euclidean space-time. 

Consider now a meson in the cm. system and take into account that for a quark one has 
£>« =► =► ipW Df=u»zV = -ifi\ (B.5) 

and p^ = iui, pp = iu 2 , so that =>• —ui, D^p =>• uj 2 , while with the 3-momentum p, 

(1) (2) 
Pi = ~Pi = Pi 



DY'^+ipi, DY'^ipi. (B.6) 
For the factor Y r (9) one obtains 

= ^tr(r(m 1 + - ipklk)T{m 2 - io 2 ^ A - iPili)) (B.7) 

Inserting T = 1, 75, 27^75, 7 M in (B.7) one arrives at the expressions (35) for different chan- 
nels. 

APPENDIX C: CALCULATION OF THE PATH INTEGRAL FOR G T , (EQ. (8)), 

IN THE FREE QUARK CASE 

To find the solution in the general case with W a 7^ 1 and an interaction depending only on 
the relative quark- ant iquark coordinates, we separate here the relative and cm. coordinates 
for any path-integral index k as follows: 

Azi - Az 2 = rj, ; = P, 

UJl + UJ 2 

L) + = UJ 1 +UJ 2 , UJ T = , (C.l) 
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with the Jacobians: 



d 3 Az l d 3 z 2 = d 3 r] d 3 p, 
dw\duj<i dui r dujj. 



Integrating (8) for W a = 1 over Y[ k dp(k) one obtains 



(C.2) 



G r (x) = 4N c Y r [ TT - ^ffl 



(C.3) 



Here the notation (D 3 f]) a i> means the initial a and final b value of relative coordinate; in our 
case, evidently, for considered current correlator a = b = 0. Also in Eq. (C.3) we use the 

1/2 



notation L 



2-irAt 



. The quantity T\ is defined as 
_ s -^,u + (k)At m 2 At u x (k)rf(k) 



2w r (fc) 

Integration over duj + {k) can be easily done, 

du + (k) exp{-£cu + (A;)f } 



n/ 



2a/cj+ - 4w r (A;) 



v / 2vrAt 



J2 k 2 Wr (fc)At 



Thus at zero cm. momentum one finds 



J G r (x)d 



3 X= I ^(^TT^#)_ e -f 2 



where 



^ = 2^(*)At + 5 ^At+!^"(*) 



2At 



At this point one can use the general relation [30, 36] 



f (D s rj 
J ~ | 



with 



2Ur % Of] 



Then for the free case, V(rj) = 0, the r.h.s. of (C.8) yields 

,2 x \ r d 3 p 







cxp 



P 

2u) r 







(2ny 



■ cxp 



P_ 

2uj r 



(C.4) 



(C.5) 



(C.6) 



(C.7) 



(C.8) 



(C.9) 



(CIO) 
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Integrating first over du T (k) and performing the steapest descent method (stationary point 
analysis) one obtains 

I 2du r (k) c ^ {k) ^iM^ = e _ 2 ^ A! (cn) 
J ^Jk).^ 



and 



/ G r (x)d*x = [ d Pj^ e -^, (C.12) 

For equal current masses m\ = m 2 = m one has evidently Q\ = lj 2 = 2u T , and from (C.ll) 
it follows that for the stationary point (for any k) 2uj t = a/p 2 + m 2 . 

To compare (C.12) with the standard Feynman amplitude for the free quark loop, one 
can go to the momentum space, 



G r (q = 0,q 4 ) = j dTe iqiT I G T {x)d 
J —00 J 

d 3 p N c Y r 



3 x 



+ 



(2tt) 3 p 2 + m 2 \2^/p 2 + m 2 - iq A 2 y /p T Tm 2 + iq 4 
d 3 p 



l.V, Vi - _ 

(2tt)Vp 2 + ^ 2 [4(p 2 + m 2 ) + q 2 } 

(C.13) 



(2tt)V(p-?) 2 

This result proves our sequence of equations for the free case. For the case of interacting 
quarks, one should use in (C.8) the spectral decomposition in the infinite set of bound states, 
as it is done in the main text, Eqs. (23,25,28). 

APPENDIX D: THE RELATIVISTIC STRING HAMILTONIAN 

We start with the relativistic SH for a meson 9192, taken in the most general form 
[22]: 

*- = E(T + ir) + r- + v ~v + f 1 w{'¥ + '*\ + £-i- L —Y < m > 

f-f 2 V 2 2uiJ 2u r J \ 2v 2) 2r 2 g(u 1 ,u 2 ,crr) 

Here u T = iO\io 2 j{uj\ + ^2) and mi{i = 1, 2) are the pole masses of a quark (antiquark) and 

g{u u co 2 , or) =u x {\- C) 2 + uj( 2 + f d(3{(3 - () 2 vd(3, (D.2) 

Jo 
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with 

C= " + > ■ (D.3) 
uji + uj 2 + or 

In the SH (D.l) taken from [22] we have added the GE potential, V GE = -f This 
can be done due to the property of additivity of the static potential in QCD [39, 40]. 
depends on the variables u>i, u 2 , and u, which have been shown to be the canonical variables 
of the SH [26], and therefore they can be defined from the extremum conditions. But first, 
instead of the operator p 2 . in (D.l) we introduce p 2 = p 2 + L 2 /r 2 and present H u as the sum 
of two terms: 

H u = H® + H BtI , (D.4) 



where 



i=l,2 

and the "string" part of the SH is 



H, 



str 



2r 2 u r 



(D.6) 



g(wi,w 2 ,r)_ 

This term occurs only for the states with L ^ and can be considered as a perturbation, 
because it gives corrections < 5% to the e.v. of the unperturbed Hamiltonian (D.5). 
(For light mesons this correction is becoming rather large only for the states with L > 5). 
Then to determine the variables u>i, u>2, and v we use the following extremum conditions 
applied to If®: 



0, ( ? = 1,2); ^- = 0. (D.7) 



duii dv 
From (D.7) it follows that 

J 2 = p 2 + m 2 , v = or, (D.8) 
and therefore can be rewritten as 



H - ] = E V m * + P 2 + ar + ^geW =T R + V (r), (D.9) 



i=l,2 



where Vo(r) is just the same potential as in (14). From (D.9) one can see that the kinetic 
term T R coincides with the one in the spinless Salpeter equation (SSE). The equation 

(iZj°> + V (r))<p nL (r) = M oVnL (r) (D.10) 
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defines the e.v. M (nL) and e.f. (p n L(r). There are two essential differences between our SH 
with the kinetic part 



T R = ^mf+p 2 + ^m 2 + p 2 (D.ll) 

and many other papers where the SSE is used. 

The first one is that in (D.ll) the quark (antiquark) mass is the pole (current) mass 
and it is not considered to be a fitting parameter. In Table I we compare the input masses 
mi,m 2 used in (D.ll) and the constituent masses from [2, 4, 5]. 

The second difference refers to the string correction H stI which can be rewritten as 

L 2 a ( gi(ui,u 2 ,r) 1 (ui - uj 2 ) 2 , C\ mi „> 



Here 



91 = Q - C + (^i + ^2 + vrfT 1 (r,u 1 ,u} 2 ), 

92 = J 1 " 1 ^,^!,^) (D.13) 



1 \ 2 ( 1 x2 



with 

J r (r 1 ,u 1 , u 2 ) = Ui {^>2 + ^crrj + u 2 + ^ar 

+ar (wi + u 2 + ar f Q - C + C 2 ) • (D.14) 

The string correction to a meson mass, A str (nL) = (H str ), can be calculated with the use of 
the expressions (D.12-D.14), in which ui, oj 2 can be replaced by their averaged values with 
a good accuracy. 

For light mesons with m-y = m 2 = 0, iO\ = uj 2 = cu, a(r) = 2u, and £ — |, the second 
term in (D.12) is absent and 

32 1 3 

Wr , W )) = T ^ (9l) = _; (92 ) = _, (D.15) 

so that in this case (mi = m 2 = 0) the string correction, 

a,,,= ( /u=- l(l+i) j< r } - ( d ' i6 > 

£>Co> 

just coincides with the string correction obtained in [24] for light mesons. It is important 
that due to the negative sign of the string correction the masses of P- and D-wave heavy- 
light mesons appear to be 30 — 50 MeV smaller in our calculations than in other relativistic 
models which use the SSE equation [2]. 
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TABLE IX: The average energies Ui(nS) = ^m 2 + p 2 ^ (i = 1, 2), the reduced mass u r , and 
the excitation energy e n (oj T ) (in MeV) for the 15 and 25 heavy-light mesons. 



Meson 




D s 


B 


B s 


B c 


ui(1S) 


507 


559 


587 


639 


1662 


w 2 (15) 


1509 


1515 


4827 


4830 


4869 


u> T (lS) 


379 


408 


523 


564 


1238 


£l(cO T ) 


541 


534 


432 


406 


149 


cji(25) 


643 


692 


741 


789 


1732 


w 2 (25) 


1585 


1590 


4862 


4865 


4898 


w r (25) 


457 


482 


643 


679 


1279 


£2(^r) 


1164 


1124 


985 


959 


687 



TABLE X: The w.f. at the origin |i? n (0)| 2 (in GeV 3 ) and (p 2 ) n s ( m GeV 2 ) in einbein approximation 
for the 15 and 25 states of heavy- light mesons. 



Meson 


15 


25 




l^i(0)| 2 


(P 2 hs 


l^ 2 (0)| 2 


(P 2 )2S 


D 


0.272 


0.273 


0.266 


0.464 


D s 


0.291 


0.290 


0.284 


0.482 


B 


0.410 


0.359 


0.410 


0.599 


B s 


0.455 


0.383 


0.439 


0.624 


B c 


1.470 


0.784 


1.032 


1.023 



In Table IX we give the values of the average kinetic energies oo\ and a>2, and the excitation 
energy e(u T ), which are needed to determine the self-energy contribution (32) to M cog (nS). 

To calculate the decay constants fp and /v we need also to know the w.f. at the origin 
(f n (0) = R n (0)/y/4n and the m.e. {p 2 ) n s- Their values are given in Table X. 

As seen from Table IX the "constituent" masses a;i(2<S') and ui(1S) of the lighter quark 
gi(orgi) differ by ~ 130 MeV for the D and D s mesons and ~ 150 MeV for the B and 
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B s mesons. This difference in ui(nS) illustrates the statement that in relativistic approach 
there does not exist a universal constituent mass uJi for different nL states, but only the pole 
(current) quark masses may be considered as universal input. 

At this point we remind that the expressions (28,35) for the decay constants have been 
derived in the EA approximation and therefore to have a consistent description one needs 
to take also |i? n (0)| 2 and (p 2 ) n s as calculated in the EA. In this approximation, instead of 
the extremum conditions (D.7), different conditions may be used, [26], namely, 

dM EA dM EA _ . 

^t^ = 0, -^ = 0, D.17 

where in the EA we use the notations Cj\ and u 2 instead of uj\ and cu 2 - The mass M EA {nL) 
satisfies the equation, 

H^(p nL {r) = M EA (nL)^ nL (r), (D.18) 
where the EA Hamiltonian [26] is given by, 

and has the same interaction V (r) and the reduced mass is given by u r = (ljilj 2 ) / + ^2)- 
Then writing 

M EA {nL) =E(| + g)+ ^ K), (D.20) 
the excitation energy e{u) r ) satisfies the equation 

P 2 1 

2^" + v o(r) J <PnL{r) = e nL (u r )(p nL (r). (D.21) 

With the use of (D.20) the extremum conditions (D.17) reduce to the equations: 

de 

Qj = m ]-2Col— (< = 1,2). (D.22) 

The derivative de(u> T ) / 'duo,, in (D.22) is a very smooth function of the variable Co r (nL) and 
can be easily calculated. Equation (D.21) formally coincides with the Schrodinger equation 
but differs from the physical point of view: the masses uJi{nL){i = 1,2) and uj T (nL)), being 
the average kinetic energy of a quark (antiquark), are different for every nL state. In this 
way relativistic corrections are taken into account in Eq. (D.21)) through the increase of Cb\ 
and l)2- 
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In the EA the values of Ui appear to be very close to Ui, determined from the SSE (26) 
or (D.10). Therefore one can define Ui(nL), solving the SSE once and calculating the m.e. 



(nL) = (vV + m 2 ) nL , (D.23) 



instead of calculations of a variety of m.e. in the solutions of (D.21) for different nL states. 
However, even for Ui = uoi the e.v. M EA (nL) (D.18) slightly differs from M {nL) for the 
SSE with the same quark pole mass m q . The differences between them may be called the 
relativistic correction 5r: 

M (nL) = M EA (nL) - 5 R , (D.24) 
which can be approximately calculated, 

5 R ~ (P 2 )- + ^--l (D . 25) 
2lo± v ; 

In (D.25) the index "1" refers to the lighter quark (antiquark), its pole mass m\ and kinetic 
energy uj 1 (D.23). Such a difference between M and M EA occurs due to the fact that 

yV + ^nL ± (sJp 2 + m 2 ) nL . (D.26) 

Surprisingly, for the D and D s and B and B s mesons 5r remains almost constant, being 
approximately equal to 70 MeV for the ground states and about 100 MeV for the 2S states. 
For the heavy B c meson 5r ~ 20 MeV is essentially smaller. 



APPENDIX E: THE SELF-ENERGY CORRECTION TO THE MESON MASS 

The self-energy correction to the meson mass originates from the NP contribution to the 
squared quark mass m 2 q as a result of the spin interaction of a quark with NP background 
gluonic field [35]. It was shown that only due to the presence of this correction in the meson 
mass it is possible to obtain a linear Regge trajectory for light mesons [24]. With the use of 
the old result from [35] and the recent result from [54] it can be presented as 

A SE („L) = £ ( -— I + + + ) . (E.1) 



i=l,2 



Here rrii(i = 1,2) is the pole mass of i-th quark (antiquark), u>i is determined by (D.7), u T 
is the reduced mass, e(uo T ) is defined by the solution of the equation (D.21). The factor rff 
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depends on the flavor of a given quark (antiquark) qi(q~i). The first, negative term in (E.l) 
was calculated in the Simonov's paper [35] while the second, positive term is rather small 
and has been derived recently [54]. The analytical expression for rjf differs for a heavy quark 
with mass m q > T g x and m q < T^ 1 , where T g is the gluonic (vacuum) correlation length 
which defines the behaviour of the bilocal vacuum correlators D(x) and D\(x) [21]. The 
value T g has been measured in lattice QCD and in the quenched approximation T g < 0.2fm 
[44] while in the unquenched case T g is larger (T g ps 0.30 fm) [45]. 

Introducing the variable y = m q T g (for c, b quarks y > 1 and for u(d), s quarks y < 1) rjf 
is given by the following expressions [35]: 

1 + 2y 2 3y 2 



arctan \/ y 2 — 1 for y > 1, 



{y 2 - l) 2 {y 2 - l) 5 /2 

Vf = T^¥?~T^¥W v ° ry< ( } 

For m u ,m d — > the quantity 77/ tends to 1, while for m g — > 00 the value 77/ — * 0. Here for 
Tfiu = iTid ~ and m& = 4.78 GeV we have 

rjf = 1.0, for u(d) quarks 

rj b = 0.03, for b quarks. (E.3) 

However, since the value of T g and the mass of the s quark are not known with a good 
accuracy, the value of rj s may essentially differ. Then from (E.2) we obtain 

T) s = 0.87 for y = 0.20 (T g = 0.22 fm, m s = 0.18 GeV), 

r] s = 0.76 for y = 0.33 (T g = 0.32 fm, m s = 0.22 GeV). (E.4) 

and also for c quarks 

r] c = 0.36 (T g = 0.2 fm, m c = 1.4 GeV). (E.5) 

In our analysis here we use the following numbers: 

Vu(d) = 1.0, Vs = 0.65, Tfc = 0.35, Vb = 0.025. (E.6) 

The contribution of the b quark to A SE (E.l)) is about —1 MeV and can be neglected. 

Finally, we notice that in the first term in (E.l) we use the number 1.5 instead of the 
number 2.0 derived in [35] (where rj u ( d ) = 0.9 was taken, as well as in [24]). This change is 
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made because here we do not neglect the contribution of the correlator Di(x) as in the first 
paper of [35]. Instead at x = we use the relation D(0) + -Di(O) = f^C 2 , where G 2 is the 
gluonic condensate (the details are given in second paper of [54]). 
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